(19) 



J 



(12) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(11) 



EP 1 429 157 A1 



EUROPEAN PATENT APPLICATION 



(43) Date of publication: 
r 16.06.2004 Bulletin 2004/25 

t 

(21) Application number: 02293082.0 
v (22) Date of filing: 13.12.2002 



(51) mt ci 7: G01V3/20, G01V3/38 



(84) Designated Contracting States: 


(72) Inventors: 


AT BE BG CH CY CZ DE DK EE ES Fl FR GB GR 


• Mathieu, Gllles 


IE IT LI LU MC NL PT SE SI SK TR 


92190 Meudon(FR) 


Designated Extension States: 


• Delhomme, Jean-Pierre 


AL LT LV MK RO 


92100 Boulogne Blllancourt (FR) 




• Cheung, Philip 


(71) Applicants: 


78110 Le Veslnet(FR) 


• Service Petrollers Schlumberger 


• Wang, Ylnyu 


75007 Paris (FR) 


78280 Guyancourt (FR) 


• Schlumberger Technology B.V. 




2514 JG Den Haag(NL) 


(74) Representative: Hyden, Martin 


• SCHLUMBERGER HOLDINGS LIMITED 


Etudes et Productions Schlumberger, 


Road Town, Tortola (VG) 


Intellectual Property Department, 




1, rue Henri Becquerel 




92142 Clamart Cedex (FR) 



(54) Method and apparatus for improved depth matching of borehole images or core images 



in 
o> 

CM 



Q. 

LU 



(57) A method for depth matching borehole images 
and / or core section images is disclosed wherein sig- 
nals from sensors at different levels on a logging tool 
are converted into an averaged signal representing the 
average bed signal at the center of the borehole at each 
of the different levels. A depth matching technique is ap- 
plied to the averaged signals from the sensors at differ- 
ent levels on the logging tool to determine the optimum 
depth offset necessary for matching two sets of signals 
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from sensors at the different levels of the logging tool. 
In an alternative embodiment of the invention a Hough 
transform is utilized to process the well log images and 
generate three-dimensional images in Hough space. 
The three dimensional images are converted into two- 
dimensional extremum curves. Depth matching is per- 
formed on the two dimensional extremum curves to cal- 
culate an offset to match the two dimensional extremum 
curves. The calculated offset is then applied to the well 
log images to depth match them. 
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Description 

Field of the Invention 

[0001] This invention relates to logging of boreholes, s 
More specifically, the invention relates to a method for 
improving correlation of core images, borehole images 
and well logging data derived from sensors located on 
a well logging tool as it passes through a borehole. 

Background of the Invention 

[0002] In a typical prior art oil well logging scenario, a 
string of well logging tools having multiple sensors for 
measuring characteristics of the earth formation along 
the wall of a borehole is lowered via a cable to the bottom 
of the borehole. Geophysical data is recorded by way of 
the sensors as the cable is wound in using a precision 
winch. The depth in the earth at which the sensors on 
the logging tools are positioned as data is logged is de- 
termined by measuring the logging speed and cable 
depth. Devices such as a depth wheel measurement in- 
strument and an axial accelerometer may be utilized. 
[0003] Typically, measurements taken along the 
length of a borehole by logging tools are intended to pro- 
vide indications of oil-bearing or gas-bearing strata in 
the earth. In the prior art, measurements of various char- 
acteristics or parameters of earth formations are usually 
obtained by combining measurements (data logs) taken 
by way of multiple sensors during a single pass through 
the borehole, or taken during different passes through 
the same borehole. When combining such measure- 
ments it is necessary that they be accurately correlated 
in depth with one another to be useful. 
[0004] One problem is that the movement downhole 
of a tool with its sensors is usually not uniform. The non- 
uniform motion may be caused by such thing as: (a) 
damped longitudinal oscillations of the logging tool on 
the cable, (b) sticking and slipping of the logging tool 
against the sides of a borehole or the wall of casing in 
the borehole, and (c) irregular motion of the sensors that 
are mounted on mechanical arms that have independ- 
ent motion with respect to the logging tool itself. As a 
result of non-uniform motion of the logging tool, data col- 
lected by a particular sensor on the logging tool at any 
specific depth in the borehole may not show as being 
recorded at the specific depth but at a different depth. 
Similarly, data from different sensors of the same tool 
may not show as having been recorded at the same 
depth for all those sensors. 

[0005] To overcome the depth measurement problem 
in the prior art, a well logging tool for measuring depth 
in a borehole is disclosed in U.S. Patent 5,019,978, is- 
sued May 28, 1 991 to Allen Q. Howard, Jr. and David J, 
Rossi. The invention in that patent provides for estima- 
tion of a dominant mechanical resonant frequency pa- 
rameter and of a damping constant parameter. These 
two parameters are taken into consideration when cor- 



recting an approximate indication of depth of a well log- 
ging tool to determine the actual, true depth of the well 
logging tool in a borehole. 

[0006] Besides measuring depth in a borehoie, well 
logging seeks to measure other physical characteristics 
along the sides of the borehole such as fractures, bed 
boundaries and bed dips. A major advance in borehole 
logging has been the development by Schlumberger of 
the Formation Microscanner ("FMS"), a borehole imag- 
ing system. The FMS processing technique is described 
in U.S. Patent Numbers 4,468,623 issued Aug. 28, 1984 
to Stanley C. Gianzero, David E. Palaith, and David S. 
K. Chan; and In U.S. Pat. No. 4,567,759 issued Feb. 5, 
1986 to Michael P. Ekstrom and David S. K. Chan. The 
FMS system uses a tool having wall-engaging pads 
each carrying an array of electrical sensors distributed 
in the circumferential direction with respect to the axis 
of the borehole. Signal voltages generated by the sen- 
sors are sampled as the well logging tool moves along 
the borehole. The signals are processed and rendered 
visible, by photographic or other printout, or by cathode- 
ray tube display as a two-dimensional visible image is 
formed over logged segments of the borehole walls. In 
such images, bed boundaries can be visually identified 
from sharp visible contrasts in the images, which reflect 
sharp changes in resistivity at boundaries of the beds. 
The images thus obtained may exhibit a resolution on 
the order of 0.5 cm, allowing very fine details of the for- 
mation to be distinguished due to the number of sensors 
in the circumferential direction , and the high rate of sam- 
pling in the longitudinal direction. 
[0007] U.S. Patent No. 4,251,773, issued Feb. 17, 
1981 to Michel Cailliau and Philippe Vincent teaches the 
use of signals from sensors to determine dip (inclination) 
and azimuth (strike) of bed boundaries. More specifical- 
ly, this patent teaches a logging tool that has four sub- 
stantially identical pads with sensors angularly distrib- 
uted about the axis of the logging tool in a side-by-side 
relationship and adapted to engage the borehole wall at 
ninety-degree intervals. The sensors provide resistivity 
measurements of the respective sectors of the borehole 
wall engaged by the pads. As the logging tool is moved 
along the borehole wall, the sensors continuously pro- 
vide signals measuring the resistivity of the adjacent 
earth formation. Sharp variations in resistivity indicate 
boundaries between different beds in the earth forma- 
tion. The signals produced by the sensors at different 
angular positions of the pads are processed to provide 
information about the dip of bed boundaries, i.e., the ori- 
entations of the bed boundaries with respect to a terres- 
trial reference, and the azimuth of the dip. 
[0008] Another tool that is lowered into or withdrawn 
from a borehole and has multiple sensors that measure 
properties such as resistivity along segments of a bore- 
hole wall is taught in U.S. Patent 5,960,371 issued Sep- 
tember 28, 1 999 to Naoki Saito, Nicholas N. Bennett and 
Robert Burridge. This patent teaches use of the Hough 
transform to extract dip and azimuth of multiple fractures 
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and beddings from any type of borehole image with re- 
spect to a terrestrial reference. The method is also ro- 
bust enough to account for noise or gaps in the images. 
The method can separate dips and azimuths of fractures 
from those of formations. Thus, it can detect and char- 
acterize other geometric features (e.g., linear, circular, 
or ellipsoidal shapes, some of which may represent 
vugs in carbonate reservoirs) present in the images. 
[0009] As previously pointed out, when combining da- 
ta of well logs it is necessary that they be correlated in 
depth in order to be useful. One method for depth cor- 
relation of well log data is taught in U.S. Pat. No. 
4,327,412 issued April 27, 1982 to John P. Timmons. 
The method disclosed in this patent is relatively complex 
and determines the displacement between a plurality of 
well logs so they may be correlated and combined. The 
well logs are derived from multiple, spaced sensors 
passed one time through a single borehole, or from sep- 
arate passes of the same sensors through the same 
borehole. First, a normalized correlation function be- 
tween selected groups of samples of the sets of data is 
determined as a first assumption to have a predeter- 
mined displacement relationship with the groups of 
samples of the data. A step of determining the normal- 
ized correlation function is repeated for a number of 
overlapping groups of well log data samples to produce 
a number of overlapping correlation functions. At least 
some of these overlapping correlation functions are 
combined to produce an improved correlation function 
that is used to depth correlate the well log data to a com- 
mon, accurate depth level. 

[0010] Another method for correlating well logging da- 
ta collected by multiple sensors is taught in U.S. Patent 
6,272,232 issued August 7, 2001 to Jean-Pierre Del- 
homme and Jean R Rivest. This patent teaches a meth- 
od for constructing, from an initial image of the wall of a 
borehole, a new "crossing-component image" centered 
on the axis of the borehole. The new image is represent- 
ative of variations in a physical parameter of the earth 
formation in both the longitudinal direction of the bore- 
hole (depth), and in the peripheral direction of the bore- 
hole wall (laterally). The new image Includes only those 
components of the physical parameters that extend all 
the way across the initial borehole image from one side 
of the image to the other. The method also includes de- 
termining variations in one or more attributes relating to 
the new image as a function of depth. The variations 
provide information relating to morphology to indicate 
solid zones, bedded zones, or different types of hetero- 
geneous zones. 

[0011] When individual well logs are concurrently ob- 
tained from a first and a second set of sensors vertically 
spaced from each other in a well logging tool the well 
logs also need to be depth correlated. They need to be 
depth correlated so that the data for any given level in 
the well taken by the two sets of vertically spaced sen- 
sors are aligned in order to be useful. Another method 
for depth correlating well log data is taught in U.S. Patent 



4,320,469, issued March 16, 1982 to William J. Frawley 
and Philip A. Mongelluzzo. This patent teaches doing 
this by using data correlograms obtained by applying a 
correlation function to a pair of digitized well logs to re- 

s duce the amount of data that must be processed. The 
remaining data is then processed in an efficient and ac- 
curate manner in order to arrive at results indicating ex- 
actly how much two logs must be shifted with respect to 
each other for an optimized depth correlation between 

10 them. 

[0012] While methods for depth correlating data ob- 
tained using sensors spaced vertically from each other 
on a well logging tool, or obtained at separate times of 
the same borehole are known, some of these methods 
15 are very complex and there is a need for a better method 
to depth correlate data obtained from well logging tool 
sensors. 

Description of the Drawings 

20 

[0013] The invention will be better understood upon 
reading the following Detailed Description in conjunction 
with the drawings in which: 

25 Figure 1 shows a well logging tool in a borehole 
along with apparatus at the surface of the earth for 
controlling the tool and recording measurements 
derived therefrom; 

30 Figure 1A is a highly simplified representation of a 
portion of a well logging tool showing the orientation 
of sensors thereon; 

Figure 2 is a picture of a portion of a well logging 
35 tool with its arms, pads and sensors; 

Figure 3A shows a bedding plane or fault intersect- 
ing a borehole at an angle relative to the axis of the 
borehole to create a dipping event and the intersec- 
40 tion is represented by an ellipse; 

Figure 3B shows a representation showing a por- 
tion of the wall of a borehole represented in Figure 
3A unrolled to a flat position; 

45 

Figure 4 shows a display of signals obtained from 
sensors mounted on two sets of pads angularly dis- 
tributed about the axis of a well logging tool, with 
each set of pads being spaced vertically with rela- 
50 tion to each other; 

Figure 5 shows signals obtained from two sets of 
sensors spaced vertically with relation to each other 
after depth correlating the signals for vertical offset; 

55 

'Figure 6 shows a block diagram of the steps in- 
volved to perform depth correlation by-a^eraging In- 
dividual ? sensor signals obtained at each level of a 
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well logging tool into a single signal, and the aver- 
aged signals are then depth correlated to determine 
the optimum depth offset necessary for precisely 
matching signals; 

Figure 7 shows averaged two dimensional signals 
obtained from two sets of sensors spaced vertically 
with relation to each other after being translated into 
three dimensional Hough space to be used to de- 
termine depth correlation to be applied to the sig- 
nals; and 

Figure 8 shows a block diagram of the steps in- 
volved in translating signals taken from the two sets 
of sensors into three dimensional Hough space to 
^ perform depth correlation. 

Detailed Description 

[0014] The terms "depth match" and "depth matching" 
used hereinafter mean correlating sets of data taken by 
different vertically spaced sensors belonging to different 
levels of sensors, such as sensors S1a - S1d and S2a 
- S2d in Figure 1 A, during a single pass of a logging tool 
through a borehole, or during different passes through 
a borehole. For a geological event, the vertical offset be- 
tween different levels of sensors causes signals coming 
from sensors in a first level of a logging tool to be gen- 
erated at a slightly different time than the signals coming 
from sensors in a second level. The difference in these 
times is also dependent on logging speed variations. 
The effects of the time-offset need to be compensated 
to "depth match" the signals coming from the sensors 
at both levels. 

[0015] Briefly, in a first embodiment of the present in- 
vention, individual traces from the sensors of any first 
specific level of sensors on a logging tool are converted 
into one signal representing those traces. The one sig- 
nal represents all the traces (averaged log signal) and 
includes an average borehole-crossing signal. The av- 
eraged log signal is generated by performing a prior art 
computation of bedding dips or flowlines (angle and az- 
imuth) on the traces from each sensor and an average 
or median computation is applied to all the samples from 
sub-sensors that comprise each sensor aligned along 
the pre-determined dip. This is repeated for a second 
specific level of sensors vertically offset from the first 
specific level of sensors and a second averaged signal 
is obtained. Then a prior art single log trace depth 
matching technique is applied to two averaged signals 
to determine the optimum depth "offset" necessary to 
precisely match the two averaged signals. The "offset" 
determined is then applied to the original traces from the 
sensors at one level to depth match them to the traces 
from the sensors at the other level. 
[001 6] Briefly, in a second embodiment of the present 
invention the Hough transform is utilized to generate 
three-dimensional images from well log signals. The 



three dimensional images are used to create two-di- 
mensional extremum curves that are in turn used to cal- 
culate an offset. The offset is used to depth match the 
well log signals obtained from vertically spaced sensors 

5 on a well logging tool, or signals obtained during differ- 
ent logging passes through the same borehole. The det- 
rimental effects of partial and non-overlapping sensor 
coverage of a borehole wall are reduced, if not eliminat- 
ed, because of the characteristics of the Hough image 

10 transform, and because depth matching is finally ap- 
plied to the two dimensional extremum curves derived 
from the three-dimensional images. More specifically, a 
depth matching offset is calculated directly from two di- 
mensional extremum curves derived from the Hough 

15 three-dimensional images without the need for calculat- 
ing dip and azimuth of bed boundaries in an earth for- 
mation and compensating for it before the depth match- 
ing may be performed. 

[0017] While the following description is given with 
20 reference to depth matching data obtained from sensors 
on a borehole logging tool traveling through a borehole, 
the teaching of the present invention may be utilized with 
records, no matter how obtained or made. The records 
may be a large number of azimuthally oriented traces 
25 obtained either from sensors applied to a borehole wall, 
from sensors scanning a borehole azimuthally, from 
sensors scanning azimuthally the cylindrical surface of 
a whole core, or from sensors probing an azimuthally 
oriented core section cut, by stacking the pixels/meas- 
30 ures along the apparent dip directions seen on the core 
section cut. 

[001 8] In the following description only two sets of da- 
ta are obtained from two sets of spaced sensors, and 
the data is correlated. This is done for the sake of clarity 

35 in describing the invention. However, the teaching of the 
invention may be used to correlate or depth match mul- 
tiple sets of data. In addition, while the sets of data col- 
lected are described as being two-dimensional, some 
types of sensors may collect three-dimensional data. 

40 [0019] In the following description of a well logging 
system, prior art material is not given in great detail to 
avoid detracting from an understanding of the two em- 
bodiments of the invention presented herein. For a more 
detailed understanding of any referenced prior art see 

45 the above Background of the Invention and the refer- 
ences cited therein. 

[0020] Referring to Figure 1 , there is shown an inves- 
tigating tool 1 0 in a borehole 1 1 for investigating subsur- 
face earth formations. Investigating tool 10 is supported 

50 in borehole 1 1 on an end of a cable 1 2. Cable 1 2 passes 
over a sheave wheel 13 and is secured to a drum and 
winch mechanism 14. The drum and winch mechanism 
14 includes a suitable brush and slip ring arrangement 
for providing electrical connections between the cable 

55 conductors and a control panel 15. The control panel 1 5 
operates to supply power and control signals to Investi- 
gating tool 10, and includes suitable electronic circuitry 
for receiving well logging signals from sensors on the 



9/18/05, EAST Version: 2.0.1.4 



7 



EP1 429157 A1 



8 



investigating tool. The well logging signals may be 
stored on a digital tape recorder 1 6. 
[0021] Tape recorder 16 converts analog signals re- 
ceived from control panel 15 into digital signals and is 
stepped as a function of depth by a driving wheel 1 7 that 5 
engages cable 12. The digital signals are then either 
transmitted or carried to a computer system 18 for 
processing. Computer system 18 comprises a proces- 
sor, display and plotter for processing and displaying 
well logging data. 

[0022] Investigating tool 1 0 has a reference point 20 
that constitutes the center or recording point of tool 10. 
Because of long elastic cable 12, investigating tool 10 
is subject to displacement arising from cable stretch 
thus causing the true depth level to be different from the 
depth level registered by a depth recorder 1 9 derived 
from driving wheel 17. The depth level registered by 
depth recorder 1 9, which is identical with that registered 
by the tape recorder 16, is designated "Z". When at- 
tempting to depth match well log data produced by dif- 
ferent sensors of investigating tool 1 0 there is a distinct 
possibility that the well logs recorded at the same depth 
level will not be referenced to the same depth level. 
[0023] On investigating tool 1 0 sensors are located on 
two vertically spaced levels L1 and L2 as shown, and 
their positions are described in greater detail starting in 
the next paragraph. These sensors produce the well log- 
ging signals that are received by control panel 15. 
[0024] Figure 1A is a simplified representation of a 
portion of tool 10 about reference point 20 showing the 
layout of sensors thereon. Sensors are usually mounted 
on pads that engage the wall of a borehole or casing in 
a borehole, although the pads are not shown in Figure 
1 A, In the highly simplified arrangement described here- 
in there is only one sensor on each pad. In Figure 1 A 
the sensors are represented by boxes with an V therein 
and are designated S1 a - S1 d in level L1 and S2a - S2d 
in level L2. In Figure 2 there is shown a more detailed 
representation of an investigating tool 10. 
[0025] Figure 1 A shows two levels L1 and L2 of sen- 
sors positioned in two parallel planes. Level L1 has four 
sensors S1 a- S1 d spaced ninety degrees apart, and lev- 
el L2 has four sensors S2a - S2d spaced ninety degrees 
apart. In one embodiment, the level L2 sensors are ro- 
tated or offset forty-five degrees with respect to the level 
L1 sensors. The spacing of the sensors in each of levels 
L1 and L2, and the rotational offset between the two 
planes, results in the orientation of all the sensors 
around a borehole as represented on the circle at the 
top of Figure 1A. 

[0026] An example of a sensor that may be used on 
tool 10 is one that measures resistivity or conductivity 
to detect and measure boundary layers in a borehole. 
Such a sensor is taught in U.S. Patent 4,786,874 issued 
November 22, 1988 to Donald S. Grosso and Allen 
Duckworth. The patented device utilizes electrical con- 
ductivity or resistivity to detect and measure boundary 
layers. 



[0027] Figure 2 is a picture of an actual well logging 
tool 1 0 with its arms 52a - 52d, pads 53a - 53d, and sen- 
sors S1 a - S1 d thereon. Arms 52a - 52d move independ- 
ently as tool 10 moves through a borehole to account 
for irregularities in the wall of the borehole. These arms 
also move perpendicularly to the axis of a borehole to 
keep sensors on the pads in contact with the wall of the 
borehole. Pads 53a - 53d are respectively connected to 
arms 52a - 52d and move therewith. On the surface of 
pads 53a - 53d are mounted sensors S1a - S1d. Each 
of the sensors S1 a - S1 d may include a number of small- 
er sub-sensors. 

[0028] When a bedding plane or fault intersects a 
borehole 11 at an angle relative to the axis of the bore- 
hole, it is Inclined or dipping and creates a "dipping 
event" and the intersection is represented by an ellipse 
57 as shown in Fig. 3a. In Fig. 3 A are representatively 
shown four traces 21a - 21 d taken by the four sensors 
S1 a - S1 d in level L1 of Fig. 1 A as the tool moves through 
the borehole. The traces are not shown extending lon- 
gitudinally along the periphery of borehole 11 to avoid 
cluttering up Fig. 3A. The four traces intersect ellipse 57 
at four points Z 1 - Z 4 . These four points represent the 
different times at which the four sensors S1 a - S1 d will 
generate a signal representing the same dipping event 
as the sensors in level L1 are passed through borehole 
11. 

[0029] When the portion of the wall of borehole 11 rep- 
resented in Figure 3A is unrolled to a flat position, the 
ellipse 57 appears as a sinusoid 58 as shown in Fig. 3B. 
The location of trace-intersect points Z A - Z 4 are shown 
on sinusoid 58. In addition, if traces 22a - 22d in level 
L2 of Fig. 1 B were shown plotted on Fig. 3A, which they 
are not since they are not in level L1 , they would appear 
in Fig. 3B as the four intersect points Z 5 - Zg and are 
interspaced between points Z 1 - Z 4 on sinusoid 58 as 
shown. With all eight intersect points shown in Fig. 3B 
it can be seen how a dipping event is represented as a 
sinusoid 58. 

[0030] When the images of a dipping event are pro- 
jected onto a flat surface, as shown in Fig. 3B, the ver- 
tical displacement of points Zj (where i = 1 through 8) 
with respect to a point Z 0 will satisfy the mathematical 
relationship: 

5Z 0 = r (<p,Z 0 + 5Z 0 ) tan (a) sin (q> -p) 

where: 

Zq = depth of the center of the event along the lon- 
gitudinal axis of the borehole; 
9 = angle between north and the imaged point; 
6Zq = Z = vertical distance between the center of 
the pattern at level Zq and a point Z; 
oc,p = dip and strike angle of the dipping event; and 
r(<p, Z 0 +8Zo) = radius of the borehole at point (<p,Z 0 
+ 5Zo)- 
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[0031] Thus, the vertical displacement distance SZg 
from point Zq to each of points Z 1 - Z 4 in level L1 may 
be calculated. 

[0032] r((p, Zq+SZq) for round cylindrical boreholes is 
constant and the borehole intersection is an ellipse. 
When the ellipse is unrolled, as described above, the 
pattern seen is a sinusoid. The amplitude of this sinusoid 
is a function of the dip angle, and its phase is a function 
of the orientation of the surface with respect to north 
(strike angle or azimuth). For a non-cylindrical borehole, 
the intersection is not an ellipse and the pattern Is a dis- 
torted sinusoid. Obviously, the amount of distortion, de- 
pends upon the degree of departure of the borehole 
from a perfect cylinder. 

[0033] Fig. 4 shows displays 21 and 22 with traces 
21 a - 21 d and 22a - 22d respectively. The set of traces 
21 a - 21 d is from the four sensors S1 a - S 1 d in level L1 , 
and the set of traces 22a - 22d is from sensors S2a - 
S2d in level L2 as seen in Figure 1A. The traces are 
obtained during logging of borehole 11. Each sensor 
may include at least one row of smaller sub-sensors. In 
addition, each sub-sensor generates a signal, and the 
side-by-side combination of such signals from the sub- 
sensors generates the eight traces 21 a - 21 d and 22a - 
22d. 

[0034] It can be seen in display 21 that sensor S1a 
produces trace 21 a at an arbitrarily assigned azimuth of 
0 degrees. With the previously described ninety degrees 
spacing between sensors, sensor S1b produces trace 
21b at an arbitrary azimuth of 90 degrees, sensor S1c 
produces trace 21c at an arbitrary azimuth of 180 de- 
grees, and sensor S1d produces trace 21 d at an arbi- 
trary azimuth of 270 degrees. The four level L2 sensors 
S2a - S2d that produce traces 22a - 22d have the same 
azimuthal spacing with respect to one another as the 
level L1 sensors, i.e., the level L2 sensors are spaced 
apart at ninety degrees among themselves. However, 
as shown in Fig. 1A, these sensors are oriented at a 
forty-five degree offset relative to the level L1 sensors 
as explained with reference to Figure 1 A. Accordingly, 
trace 22a is at forty-five degrees, trace 22b is at one- 
hundred thirty five degrees and so on. 
[0035] As logging tool 1 0 moves along borehole 1 1 , it 
may rotate horizontally a small amount and its azimuth 
orientation may change. This creates the horizontal im- 
age skewing.of the four traces shown in each of displays 
21 and 22. The skewing is shown as the traces not being 
in a vertical, straight line. 

[0036] Superimposed on display 21 are drawn two si- 
nusoidal waveforms 26 and 27 that reflect common 
borehole crossing features shown in traces 21a - 21 d. 
"Borehole crossing features" designate rock features 
that are going across the borehole as opposed to small- 
er localized rock features that may be seen by fewer 
than ail the pads of a level of sensors. Similarly, super- 
Imposed on display 22 are two sinusoidal waveforms 28 
and 29 that reflect common borehole crossing features 
shown in traces 22a - 22d. 



[0037] As is known in the prior art, and as previously 
described with reference to Figs. 3A and 3B, sinusoidal 
waveforms 26 - 29 always reflect a "dip" in a bedding 
plane or fault that intersects borehole 11. The amplitude 
s of the waveforms indicates the degree or angle of dip, 
and the location of the peaks of the waveforms (phase) 
indicates the direction or azimuth of the dip. 
[0038] As traces 21 a - 21 d and 22a - 22d of displays 
21 and 22 are generated from sensor signals taken in 
10 the same borehole, they generally would look alike ex- 
cept for a vertical displacement 30, reflecting the spac- 
ing between the level L1 and L2 sensors, and the azi- 
muth offset between the level L1 and L2 sensors. The 
vertical displacement of the traces in Images 21 and 22 
15 reflects the fact that sensors S1a - Std are physically 
located above sensors S2a - S2d on logging tool 10 as 
shown on Fig. 1 A. In addition, some vertical displace- 
ment may be caused by other factors such as stretching 
of cable 12 as previously described. 
20 [0039] Fig. 4 also shows single borehole crossing 
traces 23 and 24 that are obtained from the individual 
traces 21 a-d and 22 a-d shown in this Figure. Each sin- 
gle trace 23 and 24 respectively represents one set of 
traces 21a - 21 d and 22a - 22d, each set being convert- 
25 ed into single traces 23 and 24 respectively. In the first 
embodiment of the invention, traces 23 and 24 are ob- 
tained by performing an average or median operation 
among each corresponding set of traces 21a - 21 d and 
22a - 22d, each trace 23 and 24 representing the aver- 
so age value of a specific earth (rock) feature crossing the 
borehole. 

[0040] Each of the average traces 23 and 24 is gen- 
erated by first performing a prior art computation of bed- 
ding dips or flowlines on the individual traces from each 

35 sensor of a specific level. As a result, sinusoidal wave- 
forms, such as 26 are obtained. Then, an average or 
median computation is applied to all the sample images 
(responses from sub-sensors) located along the sinu- 
soid waveform 26 (bedding dip) obtained. By averaging 

to the signals along bedding dips, the common part of the 
signals that is due to bedding is enhanced and the lo- 
calized part that is due to heterogeneities, such as grav- 
el, vugs, etc., is attenuated. This is desirable as the com- 
mon part may be expected to be found on all levels of 

45 sensors and thus be useable for depth matching corre- 
lation. 

[0041] The determination of a sinusoidal waveform 
relies on the dip and azimuth of each bed boundary. The 
dip and azimuth of a bed boundary as a function of depth 
50 in a borehole may be determined from traces of the 
borehole by using prior art techniques such as the 
Hough transform described in the background of the in- 
vention. 

[0042] As signals 23 and 24 are obtained, a single log 
55 trace depth matching technique is applied to these sig- 
nals to determine an optimum depth "offset" for depth 
matching signal 24 to signal 23. The "offset 11 determined 
is then applied to traces 22a - 22d to depth match them 
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to traces 21a - 21 d as seen In Fig. 5. 
[0043] The averaged traces 23 and 24 may be deter- 
mined mathematically as follows. First, the intersection 
between a bedding plane and the borehole at a given 
depth Zq may be expressed mathematically as follows. 5 

Z = Z 0 + A(Z 0 ) cos ((p-<P 0 ( z o)) 

where the Greek symbols <p and Z are the coordinate 
system of the borehole image space, and A(Zq) and <p 0 
(Z 0 ) are the amplitude and azimuth of the bedding plane 
at the depth Z 0 . 

[0044] The average signal trace at depth Z 0 , S(Z 0 ), is 
then calculated as follows. Let l(q>, Z) be the pixel values 
of the traces comprising a borehole image, for each 
depth 

Z 0 : S(Z 0 ) = ^ (I ((p, Z 0 + A(Z 0 ) cos (cp - <p 0 (Z 0 )) / N v ; 

with N v being the number of valid pixels (non-absent) 
along a sinusoid at depth Z 0 . 

[0045] Depth values Zq may be obtained by way of 
the teaching set forth in the references cited in the back- 
ground of the invention (apparatus and methods for de- 
termining the required accurate depth measurements). 
The prior art teaches using such things as cables on a 
precision winch, a depth wheel measurement instru- 
ment, and sometimes an axial accelerometer whose 
output is integrated to help overcome the known prob- 
lems inherent in precisely measuring how deep a well 
logging tool is in a borehole. These depth measure- 
ments are shown as the 100 depth unit and 110 depth 
unit measurements on the graphs in Figs, 4 and 5, 
[0046] After the averaged signals S(Z 0 ) are calculat- 
ed, and displayed as traces 23 and 24, the teaching of 
the prior art is used to correlate the averaged waveforms 
23 and 24 to determine a depth offset 30 between them. 
One prior art technique for obtaining a measure of cor- 
relation to depth match two logs is disclosed in the U.S. 
Patent 4,312,040, issued January 19, 1982 to David H. 
Tinch, Bruce N. Carpenter and Elie S. Eliahou. Thetech- 
nique taught in this patent involves comparing two depth 
shifted well logs, and determining a correlation function 
indicating how closely the two logs match each other at 
a plurality of depth levels. 

[0047] The depth correlation teaching of the Tinch et 
al patent is succinctly summarized at col. 10, 1. 56 
through col. 11,1. 32. Specific detail of the correlation 
process taught in this patent is given at col. 8, 1. 16 
through col. 10, 1. 39. 

[0048] Briefly, referring to the summarization starting 
at col. 10,1 . 56 in the Tinch et al patent, two well logs A 
and B, such as from levels L1 and L2, to be correlated 
are initially considered to be depth matched in accord- 
ance with the depth indicia on magnetic tape on which 
the log data is recorded, such as the depth from depth 



recorder 19 on recorder 16 in Fig. 1. Beginning at an 
initial depth level Z a in a borehole, a first assumed value 
of the depth displacement is computed between the B 
log and the A log at each selected depth level from Z a 
to Z^, in increments of Ay. The A log is considered to be 
the base log and the displacement for the B log relative 
to the A log is computed. 

[0049] Concerning the depth correlation process at 
one depth level, the B log is effectively shifted one step 
A z at a time and a correlation function C K is computed 
at each such step using equation (1 ) shown in col. 4 of 
the Tinch et al patent. Once the correlation values of C K 
over the entire search interval 2AZ C are computed, the 
value of K which produced the maximum value of C K , 
designated K\ Is determined. The corresponding depth 
displacement will be (K'-K max/2 AZ 0 . Thus, for example, 
if AZ 0 is 1 inch and (K'-K max/2 ) is -6, the computed depth 
displacement for log B at the presently considered depth 
level Z will be -6 inches. However, this depth displace- 
ment is not computed unless C K (max) is equal to or 
greater than C K (lim). This prevents weak correlation 
quantities from causing depth displacements. If C K 
(max) is greater than C K (lim), the displacement SH(Z) 
is computed and entered into the history file for statisti- 
cal analysis. This computed value SH(Z) is only a first 
assumption of the depth displacement which may be 
changed by the statistical analysis. 
[0050] This statistical analysis takes the form of de- 
claring a depth displacement only when a consecutive 
number of identical depth displacements have been 
computed. By so doing, an occasional erroneous depth 
displacement caused by noise, for example, will not 
cause log B to be depth displaced. 
[0051] After a depth matching offset has been calcu- 
lated for waveforms 23 and 24, according to the method 
of the above cited Tinch et al patent, the calculated offset 
is applied to traces 22a - 22d shifting (depth translating) 
them upward as seen in Fig. 5. The shift is seen when 
comparing these traces in Fig. 5 with the same traces 
in Fig. 4. Traces 21a- 21 d are thus depth matched to 
traces 22a - 22d. 

[0052] When there are more than two sets of data to 
be depth matched the process is basically the same as 
described above for two sets of data. The difference is 
that one set of data (from the sensors at one depth) is 
used as a reference set of data and the offset therefrom 
to each of the other data sets is determined. The offset 
calculated for each of the other data sets is then applied 
to the other data sets to depth match to the reference 
data set. 

[0053] In Figure 6 is shown a block diagram of the 
steps involved in performing depth correlation of differ- 
ent sets of signal traces in accordance with the teaching 
of the first embodiment of the invention. At block 31 
borehole images represented by the traces 21a-21d 
from sejisors_Sja ~_ S1d are obtained and oriented azi- 
muthally. At block 32 borehole images represented by 
the traces 22a - 22d fronvsensors S2a - S2d are ob- 
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tained and oriented azimuthally. 
[0054] As previously described, the two sets of traces 
21a - 21 d and 22a - 22d reflect dip of bedding layers in 
the borehole and also reflect the azimuths of the dips. 
At block 35 the dip indicated by each of these sets of 
images is calculated using a prior art dip calculation 
technique such as the one applying the Hough trans- 
form described in the patent mentioned in the back- 
ground of the invention in connection with the Hough 
transform. 

[0055] At block 36 the calculated dip for traces 21a- 
21 d is used along with the traces themselves to gener- 
ate a single trace 23. Similarly, the calculated dip for the 
traces 22a - 22d is used along with the traces them- 
selves to generate a single trace 24. } 
[0056] At block 37 a prior art single log trace depth 
matching technique, such as taught in the above cited 
and described Tinch et al patent, is applied to traces 23 
and 24 to determine the optimum depth offset necessary 
for precisely matching the borehole traces 22a - 22d to 
traces 21a - 21 d. At block 38 the traces 22a - 22d are 
shifted relative to traces 21 a - 21 d to depth match them. 
[0057] In a second embodiment of the present inven- 
tion, two-dimensional well log traces 21a - 21 d and 22a 
- 22d of Fig. 4 are transformed into three dimensional 
Hough space using the Hough transform to get three 
dimensional images of which only one representative 
slice 39 and 40 are shown in Figure 7. The Hough trans- 
form is well known and used in the prior art to transform 
two-dimensional well logs into three-dimensional Hough 
space images from which dip and azimuth is deter- 
mined. The Hough transform as used herein transforms 
a sinusoidal waveform into a point in a three dimensional 
space and a plurality of sinusoidal waveforms making 
up the two-dimensional space into the three dimension- 
al Hough space images. These images are made up of 
a plurality of image slices 39 and 40. Slice 39 is one of 
the slices of the three-dimensional image produced by 
the Hough transform from well log traces 21 a - 21 d. Slice 
40 is one of the slices of the three-dimensional image 
produced by the Hough transform from well log traces 
22a - 22d. 

[0058] Depth matching is accurately and automatical- 
ly accomplished using the three-dimensional images 
created by the Hough transform from well logging traces 
21 a - 21 d and 22a - 22d. The three-dimensional Hough 
images have properties that make them suitable for 
depth matching of well logging images. Those proper- 
ties are: (1) insensitivity to gaps (missing data) in the 
original well logging images (the gaps in the original well 
logging images are the spaces between traces 21 a - 21 d 
and between 22a - 22d in Figs. 4 and 5), (2) insensitivity 
to rotation and orientation of the pads, and (3) dip in the 
well logging images does not have to be calculated be- 
fore depth matching of the well logging images is per- 
formed. 

[0059] Then, according to the second embodiment of 
the present invention, extremum curves 42 and 43, rep- 



resenting traces 21a - 21d and 22a - 22d respectively, 
are produced from the three-dimensional Hough imag- 
es. The extremum curves 42 and 43, reflect certain rock 
features in the earth, and represent peak values in the 
s three dimensional space at a given depth. The image 
signal peaks (local maxima) in the three dimensional 
Hough images 39 and 40 appear as peaks in the ex- 
tremum curves 42 and 43. 

[0060] Automatic depth matching is performed direct- 
10 |y on extremum curves 42 and 43 without the need for 
calculating dip and azimuth of bed boundaries in an 
earth formation and compensating for it before depth 
matching may be performed, as in the prior art. A prior 
art depth matching algorithm is then applied to the two- 
's dimensional extremum curves 42 and 43 derived from 
the Hough image slices, such as 39 and 40, to thereby 
determine a depth shift offset that will correlate extrem- 
um curves 42 and 43. The depth shift offset that is de- 
termined is thereafter applied to traces 22a - 22d to 
20 depth match them to traces 21 a - 21 d as seen in Fig. 5. 
[0061] Prior art examples of depth matching algo- 
rithms that can be used to correlate the extremum 
curves may be the same as those utilized by the first 
embodiment and may be found in: (1) U.S. Patent 
25 4,320,469, issued March 12, 1982 to William J. Frawley, 
and Philip A. Mongelluzzo, and (2) U.S. Patent 
4,312,040, issued January 19, 1982 to David H. Tinch, 
Bruce N. Carpenter and Elie S. Eliahou. 
[0062] The Hough transform utilized in the second 
30 embodiment of the invention is described in two papers: 
(1) J. Illingworth and J. Kittler, "A Survey of the Hough 
Transform", Computer Vision, Graphics, Image 
Processing 44 (1988), p. 87-116; and (2) V. F. Leavers, 
"Which Hough Transform?", CVGIP: Image Under- 
35 standing 58 (1993), No. 2, p. 250-264. See also U.S. 
Patent 5,960,371 issued September 28, 1999 to Naoki 
Saito, Nicholas N. Bennett and Robert Burridge cited in 
the background of the invention; and a paper by J. Ill- 
ingworth and J. Kittler, "A Survey of the Hough Trans- 
40 form" Computer Vision, Graphics and Image Process- 
ing, vol.44, (1988) pp. 87-116. 
[0063] How the Hough Transform can transform im- 
ages from a defined coordinate system to a different co- 
ordinate system is more fully explained in an article en- 
45 titled "Use of the Hough Transform to Detect Lines and 
Curves in Pictures", Duda, R. P and Hart, P. E. r ACM, 
vol. 15,no.1 pp. 11-15. The Hough transform technique 
was initially developed to detect straight lines in binary 
images. However, due to its simplicity, the technique 
so was extended to detect other types of simple analytic 
shapes like circles and ellipses and is used today to rec- 
ognize patterns even when they do not have an analytic 
description. 

[0064] The use of the Hough transform to automatl- 
55 cally extract features or geometric parameters from 
borehole Images was first taught in U.S. Patent No. 
3,069,654 issued Dec. 1 8, 1 962 to P.V.C. Hough. By us- 
ing the Hough transform to analyze "dip" events they 
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appear as Image signal peaks in Hough three-dimen- 
sional space. The use of the Hough transform to meas- 
ure dip and azimuth of geological beds is also taught in 
U.S. Pat. No. 5,162,994, issued Nov. 10, 1992 to David 
O.Torres. 

[0065] The Hough transform makes it possible to de- 
termine, from an image, the specific parameters char- 
acterizing a geometrical shape such as a straight line, 
a.circle, an ellipse, or a sinusoid curve, such as curves 
26 -29 in Fig. 4. The Hough transform then projects 
points of these shapes into three dimensional parameter 
space referred to as Hough space. More specifically, 
when the Hough Transform is used to process data 
points on a two-dimensional line or curve they are, math- 
ematically speaking, transformed into three dimensional 
parameter space where the dip and strike angles are 
dimensions within the parameter space. 
[0066] Generally, according to Hough formalism eve- 
ry sinusoid is a two-dimensional curve that may be 
mathematically expressed as follows. An abscissa <p, 
and an ordinate Z where Z can be expressed as: Z = Z 0 
+ A cos (<p - <p 0 ). When mathematically developing this 
equation for the Hough transform we get: Z = Zq + A cos 
(p 0 cos (p + A sin cp 0 sin <p. We then set X 0 = +A cos (p 0 
and Y 0 = +A sin <p 0 and the equation becomes Z=Z 0 + 
X 0 cos(p + Y 0 sin<p. In these equations A is the amplitude 
of a dip indicated by the two dimensional sinusoid in a 
well log image, q> 0 is the offset (phase) of azimuth of the 
dip, and Z 0 is the zero crossing depth of the sinusoid. 
The Greek symbols cp and Z are the coordinate system 
of the logging image space, and X 0 , Y 0 and Z 0 are the 
coordinate system of the three dimensional Hough im- 
age space. 

[0067] For example, by performing a summation 
along the sinusoidal waveform of the logging image pix- 
el measurements, a value may be obtained and associ- 
ated to a point (X 0 , Y 0 , Z 0 ) in the Hough space. The value 
obtained is high if the sinusoid associated to a point in 
the Hough space lies on a sinusoidal feature having a 
high pixel value in a well logging image. Conversely, the 
obtained value Is low if the sinusoid associated with a 
point lies on a sinusoidal feature having a low pixel value 
in a well logging image. In the case where there is no 
clear sinusoidal feature in the well logging image, a me- 
dium background measurement is obtained. 
[0068] The following algorithm is then used to gener- 
ate extermum curves from three dimensional Hough im- 
ages in order to apply a prior art depth matching algo- 
rithm: 

[0069] Let H(X 0 , Y 0 . Zq) be a 3D array of the Hough 
transform. 

1. Compute H average , the mean value of H(X 0 , Y 0 , 
Zq) one time for the whole array 

2. Centralize H(Xo, Y 0 , Zq) in the value domain: 

H c ( x o» Y o» Zq) = H(X 0 , Y 0 , Z 0 ) - H average 



3. For each depth Zq, find the maximum absolute 
value V(Zq) of H c : 

V(Z 0 ) = max xo Y0 { abs(H c (X 0> Y 0 , Z 0 )) at Zq } 

[0070] From the property of the Hough transform de- 
scribed above one can see that, at each depth Zq, where 
the curve V(Zq) is at a local maximum, there must be a 
dip event (either peak or trough) in the form of a sinu- 
soid, such as sinusoids 26 - 29 in Figure 4. Thus, the 
extremum curve V(Zq) characterizes the logging image 
and is a good candidate for applying a depth matching 
algorithm. 

[0071] In Figure 8 is shown a block diagram of the 
steps involved in translating signals taken from the two 
sets of sensors S1 a - S1 d and S2a - S2d into three di- 
mensional Hough space (Fig. 7) to determine a depth 
matching offset to be used to match traces 21a - 22d 
and 22a - 22d in accordance with the teaching of the 
second embodiment of the invention. At block 44 traces 
21 a - 21 d in level L1 are obtained, and at block 45 traces 
22a - 22d in level L2 are obtained, both using prior art 
apparatus and techniques as described in the back- 
ground of the invention. 

[0072] At block 46 traces 21a - 21 d in Level L1 are 
processed using the Hough transform to produce the 
three-dimensional Hough image of which one slice Is 
shown as image 39 in Figure 7. Also, at block 47 traces 
22a - 22d in Level L2 are processed using the Hough 
transform to produce the three-dimensional Hough im- 
age of which one slice is shown as image 40 in Figure 7. 
[0073] At block 48 the three-dimensional Hough im- 
age produced in block 46 is used to compute a two-di- 
mensional extremum curve, and at block 49 the three- 
dimensional Hough image produced in block 49 is used 
to compute a two-dimensional extremum curve. 
[0074] At block 50 a depth shift offset is calculated us- 
ing the two extremum curves. As previously described, 
extremum curves 42 and 43 are created from the three- 
dimensional Hough Images 39 and 40. An automatic 
depth matching calculation is performed directly on the 
extremum curves, without the need for calculating dip 
and azimuth of bed boundaries in an earth formation, to 
calculate a depth match offset. 
[0075] At block 51 the calculated depth match offset 
calculated from the extremum curves is then used to 
depth match traces 21 a - 21 d and 22a - 22d as shown 
in Fig. 5. Being as signal traces 21a - 21 d and signals 
traces 22a - 22d are offset by forty-five degrees as pre- 
viously described, the overlaid, depth matched traces 
now show combined data every forty-five degrees 
around the borehole as represented in Fig. 1 A. 
[0076] When there are more than two sets of data to 
be depth matched the process is basically the same as 
described for two sets of data. The difference is that one 
set of data (from the sensors at one depth) is used as a 
reference set of data and. The vertical offset therefrom 
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to each of the other plurality of data sets Is determined 
and then applied to the other data sets to depth match 
to the reference data set. As applied to the use of the 
Hough transform, a plurality of two-dimensional well log 
traces are each transformed into three-dimensional 5 
Hough space using the Hough transform to get three 
dimensional images. Extremum curves are then derived 
from each of the three-dimensional Hough images. Au- 
tomatic depth matching is then performed directly on the 
extremum curves, as previously described. Using the 
extremum curves resulting from one set of data as a ref- 
erence, a prior art depth matching algorithm is then used 
to compare the extremum curves representing the ref- 
erence data set with the extremum curves representing 
each of the plurality of other data sets to thereby deter- 
mine a depth shift offset that will correlate them to the 
reference set. The depth shift offset that is determined 
using each set of extremum curves is then applied to 
the corresponding data set to depth match it to the ref- 
erence data set. 

[0077] The above description describes two embodi- 
ments of a system wherein borehole logging data is col- 
lected by sensors on a logging tool passing through a 
borehole. The sensor signals representing the data are 
transmitted via a cable to the earth surface, as shown 
in Fig. 1, to be stored and processed. However, in an 
alternative embodiment of the invention, computer proc- 
essor and storage capability may be added to the log- 
ging tool and depth match processing of the data from 
the sensors is done in real time in-situ. The depth 
matched data is then either extracted from storage after 
the logging tool is extracted from the borehole, or the 
real time depth matched data is transmitted to the sur- 
face of the earth over cables that transm it data from sen- 
sors to above ground equipment as shown in Fig. 1 . Al- 
ternatively, the depth matched data from the processor 
at the logging tool may be transferred from the logging 
tool to the surface using remote telemetry capability, A 
transmitter with the logging tool transmits the depth 
matched data determined by the receiver downhole to 
a receiver on the surface of the earth without the need 
for a physical connection between the downhole logging 
tool and the receiver. Once the depth matched data is 
received by the receiver it may be transmitted to a loca- 
tion from the borehole for analysis or other processing. 
[0078] What has been described above also applies 
to images of borehole walls and borehole core, irrespec- 
tive of the measured parameter which forms the images 
such as, for example, rock conductivity, resistivity, mag- 
netic susceptibility, photoelectric factor, density of the 
formation, amplitude of acoustic reflections, photo- 
graphic brightness, etc. In addition, an image from one 
logging run in a borehole may be compared to an image 
from another logging run in the same borehole, and an 
image from a borehole may be compared to an image 
of a core from the borehole. In addition, images obtained 
from a cased borehole may be depth matched using the 
teaching of the present invention. Further, the present 



Invention may be used to depth match logs obtained by 
a logging tool from a borehole during coring as taught 
in U.S. Patent 6,003,620 issued December 21, 1999 to 
Mukul M. Sharma, Roger T. Bonnecaze and Bernard 
Zemel. 

[0079] For the sake of simplicity only, a simple ar- 
rangement for a logging tool is described. There are 
pads and sensors on two levels, each level having four 
pads spaced ninety degrees apart, and there being one 
sensor on each pad. The pads in one level are rotated 
forty five degrees with respect to the pads in the other 
level. There may be more than four pads in each level; 
the pads may be spaced other than equiangularly; the 
pads in one level may be rotated at angles other than 
forty-five degrees with respect to the pads in the other 
level; and there may be many more transducers on each 
pad. In addition, data or images from any type of trans- 
ducer, sensor, scanner or camera may be processed 
and depth matched using the teaching of the present 
invention in both of its embodiments. 



Claims 

A method for matching a plurality of data sets from 
boreholes or core sections, the data sets are ob- 
tained from sensors and are Indicative of earth for- 
mation, boundary, or interface of earth formations 
and of dip in the vicinity of the borehole, the method 
for depth matching being characterized In that: 

(a) the data sets are transformed into three-di- 
mensional images; 

(b) two-dimensional curves are derived from 
the three-dimensional images; and 

(c) an offset is derived from the two-dimension- 
al curves for applying to the two dimensional 
data sets to depth match them to each other 

2. A method for matching a plurality of data sets from 
boreholes or core sections, the data sets being ob- 
tained from sensor signals and are indicative of a 
boundary, or interface of earth formations and of dip 
in the vicinity of the borehole, the method for depth 
matching being characterized in that: 

for each data set of the plurality of data sets, 
individual signals making up the respective da- 
ta set are combined to create an averaged sig- 
nal; 

averaged signals, each corresponding to one 
data set, are processed to calculate an offset 
that correlates the averaged signals; and 
the calculated offset is applied to the data sets 
to depth match them to each other. 

3. The method In accordance with claims 1 or2 where- 
in the method is further characterized in that data 
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sets to be depth matched are obtained at the same 
time by sensors that are vertically spaced from each 
other longitudinally along the borehole. 

4. The method in accordance with claims 1 or 2 where- 
in the method is further characterized in that data 
sets to be depth matched are obtained at different 
times for the same borehole. 

5. The method in accordance with claims 1 or 2 where- 
in the method is further characterized in that a da- 
ta set to be depth matched is obtained from a core 
section. 

6. Themethodin accordance with claims 1 or2 where- 
in the method is further characterized in that the 
data sets being obtained from sensor signals are 
two-dimensional data sets. 

7. The method in accordance with claims 1 or 2 where- 
in the method is further characterized in that the 
data sets being obtained from sensor signals are 
three-dimensional data sets 

8. The method of claims 1 or 2 wherein each of said 
sensor signals is obtained from a sensor of a plu- 
rality of sensors. 

9. The method of claim 8 wherein each sensor in- 
cludes a plurality of sub sensors. 

10. The method of claim 9 wherein each signal includes 
a trace, the trace being a side-by-side combination 
of signals from the plurality of sub sensors. 

11. The method in accordance with claim 6 further 
characterized in that the transformation of the da- 
ta sets into three-dimensional images is accom- 
plished using the Hough transform. 

12. The method in accordance with claims 1 or 11 
wherein the method is further characterized in that 
the two dimensional curves have peaks Indicating 
dip events in the vicinity of the borehole. 

13. The method in accordance with claim 12 wherein 
the method is furthercharacterized in that the data 
sets have gaps in the data and the three-dimension- 
al images created using the Hough transform are 
immune from the gaps. 



10 



15 



20 



25 



30 



35 



40 



45 



50 



15. The method in accordance with claim 14 wherein 
the method is furthercharacterized in that the de- 
termined offset is applied to the data sets from the . 
vertically spaced sensors to depth match the data 
sets to each other. 

16. The method in accordance with claim 1 wherein the 
method is further characterized by the application 
of the Hough transform to depth derivatives of sen- 
sor signals, generated by sensors, to derive the 
two-dimensional curves. 

17. The method in accordance with claims 1 or 2 where- 
in the method is further characterized In that It is 
applicable to real time depth matching of data sets 
from sensors that are vertically spaced from each 
other longitudinally along the borehole. 

1 8. The method of claim 2 wherein said averaged signal 
obtained by determining an average of the sensor 
signals along the bedding dip for a given depth in 
the borehole. 

19. The method of claim 18 wherein said computation 
of bedding dips for the sensor signals is performed 
by way of the Hough transform. 

20. A well logging tool comprising a plurality of sensors, 
a processor, and storage means coupled therebe- 
tween, the plurality of sensors adapted to generate 
data sets indicative of earth formation, boundary, or 
interface of earth formations and of dip in the vicinity 
of the borehole, the well logging tool performing the 
method of claim 1 . 

21. A well logging tool comprising a plurality of sensors, 
a processor, and storage means coupled therebe- 
tween, the plurality of sensors adapted to generate 
data sets indicative of earth formation, boundary, or 
interface of earth formations and of dip in the vicinity 
of the borehole, the well logging tool performing the 
method of claim 2. 

22. The well logging tool of claims 21 or 22, said logging 
tool adapted to transmit the depth matched data 
sets to the surface. 



14. The method in accordance with claim 1 wherein the 
method is further characterized in that two-dimen- 
sional curves for data sets from sensors that are 
vertically spaced from each other longitudinally 
along the borehole are processed to determine an 
offset that will match the two-dimensional curves. 
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